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1.  INTRODUCTION 


This  report  is  a  compilation  of  considerations  involved  in  the  design  of 
elect ro -magnetic  rods  for  satellite  erection.  Formulae  are  derived  show- 
ing  the  effect  of  various  choices  of  core  and  wire  material  and  weight, 
number  of  cores,  length  of  cores,  number  of  turns  per  core,  and  wire  area 
on  the  magnetic  moment,  hysteresis  power,  hold  power,  hold  voltage,  and 
switching  time.  These  formulae  are  then  applied  to  the  problem  of  produc¬ 
ing  a  required  magnetic  moment  at  the  least  cost  in  weight  and  power  con¬ 
sumption. 

2.  DERIVATION  OF  CORE  AND  WINDING  RELATIONS 

The  following  sections  are  devoted  to  the  derivation  of  those  equations 
that  describe  the  functioning  of  a  single  core  with  a  single  winding.  Subse¬ 
quent  sections  will  treat  questions  of  multiple  rods  anchor  windings. 

2.  1  Determination  of  Fields  in  Core 


In  this  section  we  derive  the  dependence  of  the  core  fields  on  the 


drive  current,  etc.  This  simple  derivation  makes  certain  gross  approxi¬ 
mations,  but  is  justified  by  the  more  complete  analysis  in  APPendix  1.  The 
core  geometry  and  relevant  notation  are  introduced  in  Figure  1. 

We  begin  by  assuming  a  uniform  axial  magnetization,  M.  This 
magnetization  charges  the  rod  ends  to  a  surface  density  of  /u  We  shall 

treat  the  end  charges  as  points  with  total  charge  ^MA  and  compute  the  re¬ 
sulting  demagnetizing  H  along  the  Z  axis.  It  is  easily  verified  that 


2  MA 


1  f  (2z/  £  r 


1 


3-63-3^7 


SURFACE  CURRENT  h  =  Ni  /  J 


i - 

H  -  TOTAL  H  IN  CORE 

H  -  DEMAGNETIZING  H  DUE 
d  TO  MAGNETIC  CHARGE 

M  CORE  MAGNETIZATION 

8  CORE  B 

f  -  NUMERICAL  CONSTANT 
CV  =  A  t  ~  CORE  VOLUME 
N  TURNS/WINDING 

CURRENT  IN  WINDING 


Figure  1  Core  geometry  and  notation 
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is  clearly  non-uniform;  it  is  a  minimum  at  the  center  (z  =  0)  and  tends 
to  infinity  at  the  ends  (z  =  ±  l/Z)  of  the  rod.  It  is  convenient  to  introduce  an 
"average"  and  employ  it  as  though  the  demagnetization  field  were  uni¬ 
form.  Hence,  we  shall  supplant  the  previous  relation  by 


H 


d 


f  MA 
> 

r 


f  M  CV 

- 5 

l 


The  numerical  factor  f  must  be  determined  experimentally  or  by  a  more  in¬ 
tricate  analysis  such  as  that  in  Appendix  1.  It  generally  lies  between  three 
and  ten. 

The  applied  current  produces  a  nearly  uniform  field  in  the  core 
region  of  magnitude  k.  The  total  H  field  is 


H 


We  can  eliminate  M  from  this  expression  by  employing 


H  f  M 


This  gives  the  load-line  equation 


i  - 


f  CV 


) 


H  + 


f  CV 


k 


Since  we  are  interested  only  in  long  thin  rods 

(  cv  1 

1  “M"  <  TDM 

i 

and  so  this  ratio  may  be  ignored  in  the  coefficient  of  H. 

The  load-line  can  be  combined  with  the  B -H  curve  for  the  material 
to  obtain  the  operating  point  of  the  material.  Figure  Z  shows  this  construc¬ 
tion. 
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LOAD  LINE 
WITH  THE  EQUATION 


Figure  2  Construction  for  determining  operating 
point  at  core- 
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Z.  Z  Determination  of  the  Magnetic  Moment  of  the  Saturated  Rod 


Outside  the  core,  the  effect  of  the  current  sheet  is  identical  to  that 
of  a  volume  dipole  moment  of  density  k.  Hence  the  magnetic  moment  i'O  is 

<$V\  =  CV  (k  +  M) 

We  already  have  expressions  for  k  and  M  in  terms  of  B  and  H.  Inserting 
these  gives 


r ,  .  r  cv  ^  b  f  cv  H 


™  =  <=VL0  1  j 


l 


For  saturated,  high  permeability,  long,  thin  rods  we  approximate 


where  B  is  the  saturation  B -field  of  the  core, 
s 


Z.  3  Determination  of  the  Currents  Required  for  Saturation  and  Holding 

The  required  saturation  and  holding  currents  depend  on  the  core 
material.  Figure  3  shows  a  typical  B  -H  curve  and  introduces  some  nota¬ 
tion.  From  Figure  Z 

k 


Hence:  i 

s 


Hi  = 


=  H+  fCV  _B_ 


saturation  current  =  k?£/N 


hold  current  =  k^  jP /N 
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s 

AH  -  ' 'WIDTH"  OF  LOOP 

H,  =  H  USED  FOR  HOLDING.  H  MUST 

AT  LEAST  EQUAL  H  —  AH 

s 


Figure  3  A  typical  B-H  curve  and  some  relevant 
parameters 
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where  we  have  defined  the  two  parameters 


k,  =  f  ff'rtVi’3 
k  ^  =  Hh  f  f/VV  /3 

Z.4  Determination  of  Hysteresis  Loss  in  Core 
The  hysteresis  loss /cycle  is 

hyst.  loss/ cycle  =  ZB  AH  CV  =  Z  fj.  Of]  A  H 
'  s  o 

Hence,  if  the  length  of  a  cycle  is  T,  the  hysteresis  power  dissipation  is 

ZuJSO  A  H 

P  =  - - - 

hyst  T 

Z.  5  Determination  of  Switching  Inductance 

The  switching  inductance  is  computed  by  dividing  the  flux  linkage 
change  by  the  current  change  during  traversal  of  the  non  -  saturated  portion 
of  the  B-H  curve.  From  Figure  3,  the  change  during  switching  in  B  and  H 
are 

change  in  B  =  Z  B^ 
change  in  H  =  Z  H  -  AH 

Letting  N  equal  the  number  of  turns,  we  see 

flux  linkage  change  =  Z  AB^N  =  Z  fjJ) Y'}  N/  £ 

Since  we  have  an  expression  for  current  in  terms  of  B  and  H  we  can  com¬ 
pute  the  current  change: 
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current  change 


N 


H 


A  H  t 


f  cv 


T 


M. 


2  f 
"N" 


hi 


A  H 


t 


f 
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Hence  the  switching  inductance  L.  is 


2.b  Wire  Area  and  Coil  Resistance 

The  winding  geometry  and  notation  are  shown  in  Figure  4.  We 
shall  express  the  resistance  of  the  coil  in  terms  of  WV,  CV,  N,  f  and  o. 
F rom  the  figure 


t  WV 


7T  t.  R  -  CV 


Hence 


R 


CV  + 


V  7T  p 


F  WV  J 


1/2 


But 


G  GV  'V2 

r  '  b  TT ) 


The  mean  turn  length  is  tt  (R  f  r)  so  the  wire  area  is: 


A  -  WV 

Awire  N  tt  (R  +  r )  ~~  ,  T 


WV,)1//n(^) 


L1  f  V*  1  "cv-  J  J 


It  is  convenient  to  introduce  a  parameter  defined  by 

u  .  /  “  CV^[,  ,?WV  y 1 

ki  v-  wv  )  l  1  +  v  1  f  cv  y 
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WV  CONDUCTOR  VOLUME 

*  PACKING  FACTOR -THE  RATIO  OF  WINDING 
VOLUME  TO  CONDUCTOR  VOLUME 

N  NUMBER  TURNS 
p  =  CONDUCTOR  RESISTIVITY 
l<1  A  DEFINED  PARAMETER 


Figure  4  Winding  geometry  and  Notation 


In  terms  of  this  parameter  we  may  write 


A  . 
wi  re 


1  /2 

r wv ; s  7  i 


A  ? 


N 


The  total  length  of  the  vire  is  WV/A  .  so  the  resistance 

wire 


R  „ 
coil 


WV 

P  -5- 

A"  . 
wire 


5  kjN 

=  P— 7— 


Z.  7  Hold  Voltage  and  Power 

The  hold  voltage  and  power  are  easily  expressed 
other  parameters.  The  hold  voltage  is 


v'  =  !h  Rcoil  =  P?k,  k,N 


hold 


1  3 


The  hold  power  is 
P 


,  ,  ,  =  i ,  R  .,  =  k.  kf  l 

hold  h  coil  1  3 


For  future  use  we  note  the  identity 


A  .  V,  .  , 
wire  hold 


A 


1/2 


WV  P  Phold/ 


Z.  8  Voltage  Source  Switching  Time  Constant 

^he  time  constant  appropriate  for  switching  from 


source  is 


switch  R 


fj,  <r) 

o 

p=  ki  vz  2  ) 


\  of  the  coil  is 


in  terms  of  the 


an  ideal  voltage 
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2.9  Single  Rod  Formulae  -  Summary 


fffi  = 

CVB  /[x 

R  coil  = 

P§  kjN lZ l 

where  = 

'Nr, 

vrwvA1  ’ 

i  - 

s 

k  j?/N 

whe  re  k  = 

H  + 
s 

f(VN/^3 

k3  P./N 

where  k  ^  = 

Hh * 

^hold 

0?  k^N 

P  = 

hy  st 

2  mo<T0a  h/t 

Phold  = 

p-klk3  p 

Awire  Vhofd  (WVCPhold»1/,! 

L 

<kz  - 

A  H/2) 

T  .  = 

switch 

M0<no/p?  kj 

Mk2 

,  -  AH/2) 

.  §  WV  ^1/2  "1  2 

+  cv  J 
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3.  TRADE  OFFS  AND  OPTIMIZATIONS 


This  section  applies  the  formulae  derived  in  Section  3  to  the  problem  of 
design  optimization.  The  parameters  that  the  designer  can  vary  are: 

Magnetic  material 
Winding  material 
Core  volume 
Wire  volume 
Core  length 
Wire  area 
Number  of  turns 
Hold  voltage 


Not  all  of  these  are  independent,  and  certain  of  them  may  be  fixed  by  ex¬ 
ternal  constraints.  These  parameters  must  all  be  selected  to  produce  a 
given  (5 which  is  generally  determined  by  the  dynamics  of  the  system  . 
Moreover,  the  designer  wants: 


12 


Minimum  weight 
Minimum  hysteresis  power  loss 
Minimum  hold  power  loss 
Minimum  switching  time 

He  is  limited  by  the  practical  constraints  that: 

Wire  size  must  not  be  too  small 

Core  length  must  not  exceed  the  satellite  diameter 

3.  1  Determination  of  Core  Material 

Once  ^  has  been  determined  and  a  material  selected,  the  hyster¬ 
esis  power  loss  is  fixed  independent  of  all  other  design  choices.  The  mate¬ 
rial  with  lowest  A  H  also  has  the  lowest  hysteresis  loss.  However,  the  hys¬ 
teresis  loss  is  often  a  small  fraction  of  the  hold  power,  and  the  dependence 
of  hold  power  on  core  material  is  somewhat  more  complex.  The  choice  of 
material  affects  the  hold  power  two  ways: 

(i)  Materials  with  high  B  require  less  material  to  produce  a 

s 

given  jn .  Thus,  for  a  fixed  system  weight,  more  wire  can 
be  used,  reducing  . 

(ii)  Materials  with  a  low  H^  will  have  low  values  of  k?. 

Clearly,  if  we  could  find  a  material  that  has  highest  B  and  lowest  H.  ,  it 

s  h 

would  be  the  best  choice.  Since  so  few  materials  are  even  contenders,  the 
simplest  approach  is  to  work  through  the  design  procedure  to  be  given  later 
for  each  and  select  the  one  that  best  meets  the  mission  requirements  as  to 
weight  and  hold  power. 

3.2  Selecting  Core  Length  to  Minimize  P^old 

As  we  have  seen  in  Section  2,  the  hold  power  is  proportional  to 

^  R3  A2 

phoid  ~  1  +  tt)  ' 


1  3 


hold 


<  i _ i _ i - 1 - 1 - 1— 

0  1  * 

i/0 


4 


5 


Figure  5  Graph  showing  how  Pnold  depends  on 
core  length. 
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where  we  have  just  defined 

83  =  f  l'VHh 

Differentiating  this  with  respect  to  /  gives 


d  P 


hold 


d  p 


) 


Thus,  as  P  increases  from  zero,  P,  decreases  to  a  minimum  and  then 

hold 

rises  steadily.  This  minimum  value  obtains  when 


/  = 


S 


If  the  rod  length  is  constrained  to  be  less  than  some  P  we  should  select 

max 


/  =  ,!/3 


inff  5  '  B,  £  ) 

v  max/ 


in  order  to  minimize  the  holding  power.  Figure  5  displays  the  dependence 
of  the  holding  power  on  core  length. 


3.  3  Effect  of  Replacing  a  Single  Core  by  Several  Small  Ones 

Under  certain  circumstances  it  is  advantageous  to  replace  a  single 
rod  by  several  smaller  ones.  The  total  core  colume  must  remain  the  same 
in  order  to  produce  the  required  magnetic  moment.  Moreover,  if  the  weight 
is  fixed  the  ratio  WV/CV  for  each  core  is  unchanged.  However,  the  terms 
P  do  change,  and  a  significant  reduction  in  hold  power  may  result.  Sup¬ 
pose,  for  example,  we  choose  to  employ  e  rods.  Then  from  Section  Z,  the 
total  hold  power  is  proportional  to: 


hold  -  total 


e 

e* 


1  + 


e  P. 


max 

"3 


max 
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where  we  have  defined: 


fflfl/H.  r 

h  max 


p3//3 


max 


and  ))  is  the  total  magnetic  moment  of  the  system.  We  have  seen  in  the 
previous  section  that  the  core  length  should  be  selected  as 


intl  vr) 


1/3 


max  max  J 


«  e 
r  or  — 

£ 

>  5  we  have  f 

<  £  and 

max 

: 

p 

hold -total 

, 

~  /.ij 5 

1/3  (  eV^3 

u; 

When 

e  r 

—  <  5  we  see  &  - 

£ 

and: 

max 

p 

hold -total 

~  f  -  f  2 
V  e 

O) 

This  latter  has  a  minimum  at  —  =  1.  Figure  6  shows  how  t  le  hold  power 
depends  on  the  ratio  e/e.  The  best  choice  of  e  would  be  e  =  e  but  this 
may  not  be  an  integer.  We  see  immediately  that  splitting  is  never  profit¬ 
able  if  £  <  1 .  If  splitting  is  profitable,  it  should  never  proceed  past 
e  =  £  f  1  ^  Z  c  <  5  c.  Hence,  we  need  only  consider  the  portion  of  the 

curve  in  Figure  6  corresponding  to  £  -  l 

&  r  to  max 


We  note  that  e  t  1  rods  require  less  hold  power  than  e  if  and  only 


if: 


£  >  1 
e  <  £ 


2-Li  t  z  *  *  <  *  *  i  y  £ 

e  e  +  1  e  e 
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Figure  6 


Graph  showing  how  total  hold  power 
depends  on  number  of  cores  used  to 
produce  a  given  magnetic  moment 
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Figure  7  Graph  showing  the  dependence  of  hold 
power  on  wire  volume. 
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RELATIVE  WIRE  GAUGE 


£  (  WV  /cv) 


Figure  8  Graph  showing  how  wire  gauge  depends 

on  wire  volume  for  constant  hold  voltage, 
the  effect  of  varying  wv  is  to  change  the  true 
wire  gauge  by  the  change  in  relative  wire 
gauge  determined  from  this  graph0 
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These  are  all  equivalent  to  the  statement: 


e  (e  t  1 )  < 


Hence,  we  have  the  following  role  to  determine  the  optimum  integral  number 
of  rods: 

e  =  1  if  e  ^  /T7T 

opt 

Z  if  /ITT  <c  £  <.  T273 

3  if  /zn  <  c  ^  /ITT 

etc . 

For  any  choice  of  e  between  1  and  e  ^  the  core  lengths  are  f  -  f  .So 

opt  max 

long  as  e  is  smaller  than  5e  we  have  the  convenient  formula: 


hold -total -e  cores 
hold-single  equivalent  core 


e 

G 


f  z 


f  z 


f 


f 


c 

e 


G 


3.4  Wire  Volume  Considerations 


In  Section  Z  we  saw  that  the  hold  power  was  proportional  to  k ^ . 
Figure  7  shows  the  dependence  of  on  the  ratio  §  WV/CV.  The  asymptote 
as  WV  -*  00  is  at  7T  .  From  this  graph  we  see  that  there  is  little  point  to  us¬ 
ing  more  wire  volume  than  10  CV/^f  unless  power  is  at  an  extreme  premium. 
(With  aluminum  wire  and  a  ferromagnetic  core  this  corresponds  to  a  wire 
to  core  weight  ratio  of  approximately  Z.  4).  However,  it  is  true  that  more 
wire  always  reduces  the  hold  power. 


It  is  also  true  that  increasing  the  wire  volume  may  ease  the  wind¬ 
ing  problem  by  permitting  larger  wire  area  for  a  given  hold  voltage.  Fig¬ 
ure  8  shows  the  dependence  of  wire  area  times  hold  voltage  on  ware  volume. 
It  is  based  on  the  result,  easily  derivable  from  the  formulae  in  Section  Z, 
that: 
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A  .  V,  ,, 
wire  hold 


1  f 


1  1-5  WV  y 

CV  J 


1/- 


and  the  wire  gauge  formula: 

Gauge  Z  -  Gauge  1 


Log  ( 

V  Area  Z  / 


Log  2 

3.  5  Winding  Material  for  Minimum  Pj^ 

In  Section  Z  we  saw  that  the  hold  power  was  proportional  to: 


Phold  p  kl 


There  are  only  two  contenders  for  the  wire  material  -  copper  and  aluminum. 
Their  densities  and  resistivities  are  shown  below: 


P  (A  1 )  = 

2.83*10  8  r-m 

d  (A  1 )  = 

2.  7-  10J  kg/m  5 

P  (Cu)  = 

1 . 72*  10"8  r-m 

d  (Cu)  = 

8.  9-  10J  kg/m3 

As  is  shown  in  Figure  7,  for  large  amounts  of  wire  approaches  asympto¬ 
tically  to  7i  .  Under  these  conditions  the  better  conductor,  copper,  mini¬ 
mizes  the  hold  power.  On  the  other  hand  for  low  permissible  wire  weights, 
aluminum  minimizes  the  hold  power  since  the  increase  in  wire  cross-section 

area  more  than  compensates  for  the  increase  in  resistivity.  We  shall  assume 

3  3 

a  ferromagnetic  core  with  a  density  of  8.  13*  10  kg/m  .  We  shall  let: 


Then 


r\ 


WV 

CV 


ratio  of  wire  weight  to  core  weight 

core  density 


wire  weight 
wire  density 


core  weight 


=  *n 


8. 13- 10' 


wire  density 
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Hence,  the  ratio  of  the  hold  powers  with  aluminum  and  copper  wire*  is: 


phoid(A1)  _  i  r  (i  +  3. 01  f  f  i  n2 

Ph  )ld(Cu)  3  L  (1  f  .  914  i  n,)l/'3  f  1  J 


The  cross  over  point  where  A1  and  Cu  yield  equal  hold  powers  occurs  when: 

r-r  _  1  f  (1  +  3.  01  F  n)1//2 

1  f  (1  +  .  914  ?  t|  ) 1 7 

This  occurs  when  ^  r\  -  2.83.  Hence  for  Z  r)  <  2.83  aluminum  is  the  better 

choice.  We  shall  see  that  aluminum  also  permits  the  use  of  larger  wire. 
Thus,  it  may  be  necessary  to  use  aluminum  even  when  copper  would  give  a 
slightly  lower  hold  power.  Even  when  £  r|  =  102,  the  hold  power  for  alumi¬ 
num  exceeds  that  of  copper  by  only  50  per  cent.  In  no  case  is  aluminum 
worse  than  copper  by  more  than  65  per  cent. 

3.6  Selection  of  Wire  Material  for  Maximum  Wire  Area  at  a  Given 
Hold  Voltage 

Using  results  derived  in  Section  2  we  see  that: 


wire 


A  . 
wire 


(Al) 

(Cu) 


p  P  (Al)  Phold  (A  1  )/d  (Al)  -|l/2 
L  p  (Cu)  Ph(jld  (Cu)/ d  (Cu)  _ 


=  Z.  33 


r  phoid (A1)  yA 
L  phdd  <c"> J 


As  seen  in  the  previous  section,  the  ratio  of  the  two  hold  powers  always 
liesbetween.5for5'n  =  0  and  1 . 65  for  %  T]  =  Hence ,  the  area  ratio 

always  lies  between  1.65  and  2.99.  Thus  aluminum  always  permits  the 
use  of  larger  wire.  The  difference  in  wire  gauge  numbers  ranges  from 
2.  17  to  4.  75. 
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3.7  Voltage  Source  Switching 


One  way  to  reverse  the  current  in  a  rod  is  to  apply  a  back  voltage 
to  the  entire  coil  as  shown  in  Figure  9.  Zero  time  is  required  for  the  op¬ 
erating  point  to  move  along  the  saturated  portions  of  the  X  -  i  curve,  so 
switching  actually  begins  at  a  current  of  -  i  t  Ai.  Ai  is  the  width  of  the 
X  -  i  curve  --  if,  as  we  shall  assume,  i^  ,  is  selected  to  hold  on  the  knee, 
then  Ai  =  i  -  i^.  In  any  event,  Ai  =  %  AH/N. 

Figure  10  shows  the  switching  transient.  From  it  wre  can  compute 


V 

s 


+ 


s 


1 


e 


-t/r 


Evidently  the  switch  voltage  must  exceed  the  hold  voltage.  By  direct  inte¬ 
gration  one  may  verify  that  the  total  charge  drawn  from  the  source  is: 


Q 


which  can  be  written 


Q 


(r  f  1 


-  (1  -  e' 
(1  -  e  ~x) 


n 


i 


whe  re 

x  =  t/j 

r  =  is/ih 


For  x  ^  Z  the  quantity  in  brackets  is  closely  (within  a  few  per  cent)  ap¬ 
proximated  by  .5  f  .08x.  Hence,  for  —  ^  Z  the  charge  drawn  from  the 
switching  supply  is: 


Z3 


3-63-3855 


x 


THE  DIODE  JUST  HOLDS 
THE  INDUCTOR  CURRENT 
DURING  SWITCHOVER 


igure  9  Voltage  source  switching 
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3-63-3856 


Figure  ID  Corrent  vs  time  during  voltage  source 
switching. 
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( 


1  “  xu 

s  h 


.  5  f  .  08  -  +  .08  ! 

T  y  T 


The  energy  drawn  from  the  supply  is  of  course  V  Q .  Note  that  when  the 
switching  source  is  supplying  power,  the  hold  source  is  not  so  that  this 
drain  does  not  represent  additional  power. 

3.8  Summary  of  Design  Formulae  and  Results 


l  =  inf  (j> 


,1/3 


B, 


ma 


J 


whe  re 


f  /Hu 
h 


e 

opt 


1  if  e  ^  ym 

2  /TTT  L  e  <  ,/TTT 

3  /27T  L  e  ^  /TT4  etc. 


where  € 


3 

max 


P 

P 


hold  -e  cores 
hold-1  core 


e  /g  +  2  -f  c/  e 

i  7  c  +  2  +  e 


A1  gives  a  lower  hold  power  than  Cu  for  low  wire  weights.  The  cross-over 
occurs  when 


_  wire  weight  .. 

%  - r^T —  =  2.  83 

core  weight 

A1  always  leads  to  larger  wire  size  for  fixed  hold  voltage.  The  difference 
is  on  the  order  of  3  gauge  units. 
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T7=^ 


J 


l. 
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i  t 
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rl  s 


- 1 


-  ^ 


5  +  .  08 


t) 


f  .  08  - 
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t_ 

T 
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APPENDIX  A 


The  Wound  Prolate  Ellipsoid 


1.  A  Prolate  Ellipsoid  Uniformly  Magnetized  Along  Its  Axis  of  Rotational 
Symmetry 


In  this  section  we  shall  determine  the  field  resulting  from  an  ellipsoid 
with  uniform  magnetization  along  its  rotary  symmetry  axis.  Figure  A-l 
shows  the  geometry 


Figure  A-l  The  Prolate  Ellispoid 

The  effect  of  the  uniform  magnetization  is  to  produce  a  surface  charge 
on  the  ellipsoid  with  density 


K„  i  z  .  ,  n 

a  =  /j  Q  Ml  1 

One  can  easily  obtain  an  explicit  form  for  a  and  verify  the  following  inter¬ 
esting  fact:  if  the  charge  is  all  projected  normally  onto  the  z-axis,  the  re¬ 
sulting  line  charge  density  is  proportional  to  z  --  and  the  total  positive 
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charge  equals  jj.  times  the  area  of  the  ellipsoid  at  its  equator.  Figure  A- 
shows  this  line  charge  %. 


£  =  DENSITY  OF  CHARGE  3-63-3858 

PROJECTED  ONTO  2  AXIS 


Given  the  charge  one  can  find  the  fields  directly  by  Greens -function 
methods,  or  else  one  can  solve  for  the  magnetic  potential  by  inserting  the 
charge  as  a  boundary  condition  and  expanding  in  ellipsoidal  harmonics. 

The  following  two  facts  emerge:  first,  the  internal  H-field  is  uniform;  and 
second,  the  external  field  is  the  same  as  that  produced  by  a  line*  of  charge  , 
linearly  graded,  extending  between  the  two  foci,  with  dipole  moment  equal 
to  M  times  the  volume  of  the  ellipsoid.  Figure  A -3  show  s  this  equivalent 
density 


C  =  CHARGE  DENSITY  PRODUCING  3-63-3859 

SAME  EXTERNAL  FIELD 


Figure  A  -  3  Charge  Producing  Same  External  Field  as  Ellipsoid 
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We  can  determine  a  as  follows: 


[  z-  azdz 

O 


4  7T  a  b  M 


Performing  the  integration  and  solving  for  a  gives: 

.  x/r  a  b 

a  =  Z  7T  fJL  M  - 5 - y  -v  /  V 

°  -b-)S  - 

Now  that  a  has  been  determined  we  can  compute  the  external  H-field  at  the 
equator.  Figure  A -4  shows  the  relevant  geometry. 


IV 

dH  = 

azdz 

4  TTp0d2 

a zd  i  i  / 

X 

4  w  u  d2  d 

V 

^  <) 

b 

N\C|  =  U2+  I; 

i - - 

* 

— i - .Xl - . 

“  v/a2-  b2  ^ qZ  - 


Figure  A -4  The  Determination  of  H  at  the  Equator  of  the  Ellipsoid 


By  symmetry  it  is  evident  that  H  will  have  no  normal  or  rotary  component 
on  the  equator.  From  the  figure,  the  value  of  the  H-field  in  the  negative  z 
direction  is: 
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a 


H 


r~> - > 

a  f  b" 


4  7T  jU 

o 


3 

z 


dz 


.  /n 72 

a  t  a  f  a  -  b 

Log  -  -  ?. 


4  TT  U  V 
o 


r~i  t< 

a  -  a  -  b 


Jl  71 

a  -  b 

a  J 


Since  the  tangential  component  of  H  is  continuous  across  the  boundary,  the 
uniform  H  field  within  the  ellipsoid  points  in  the  negative  z  direction  and  has 
magnitude  equal  to  the  above.  Thus  the  internal  H  opposes  M  and  the  mag¬ 
nitudes  are  related  by: 


H  =  B  M 


whe  re 


ab 

L 


Log 


a  +/”2  72  zy~l  72-, 

a  -  b  a  -  b 


a  -f~L  72 
a  -  b 


If  we  let  A  =  a/b  we  get: 


B  = 


A  -  I 


Log  Aa  f  W  -  1 J  -  1 


3.  Fields  Due  to  Current  Circulating  About  a  Prolate  Ellipsoid 

To  find  the  fields  we  exploit  the  equivalence  of  circulating  current  and 
dipole  sheet.  The  H  field  due  to  current  filament  i  equals  (except  on  the 
sheet)  that  due  to  a  sheet  of  dipoles  with  surface  density  i  that  spans  the 
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filament.  As  shown  iu  Figure  A -5  a  surface  current  whose  projection  on 
the  z-axis  is  uniform  density  k,  produces  the  same  external  field  as  an 
ellipsoid  with  uniform  magnetization  k. 

3-63-3861 


Figure  A -3  Dipole  Sheet  Construction  for  Determination  of 
External  Field 


A  similar  approach  can  be  used  to  find  the  field  inside  the  ellipsoid  In 
this  case  the  dipole  sheets  must  not  penetrate  the  interior.  Figure  A -6 
shows  the  appropriate  construction. 


3-63-3862 


Figure  A~6  Dipole  Sheet  Construction  for  Determination  of 
Internal  Field 


32 


Thus  the  internal  field  is  that  produced  by  an  infinite  slab  of  uniformly  mag¬ 
netized  material  with  an  ellipsoid  hole  in  it.  As  shown  in  Figure  A -7,  this 
field  is  easily  determined  by  superposition 


3-63-3863 


Figure  A -7  Determination  of  Internal  Field  by  Superposition 


The  slab  produces  an  H  field  in  the  plus  z  direction,  since  it  behaves  like  a 
parallel  plate  capacitor  with  surface  charge  density  [j,  k  .  The  magnitude  of 
this  field  is  k.  The  ellipsoid  produces  a  field  of  0  k  in  the  negative  z  direc¬ 
tion.  The  total  internal  field  is  thus  (1-  3)  k  in  the  plus  z  direction. 


3.  Determination  of  Fields  in  a  Wound,  Magnetizable,  Prolate  Ellipsoid 
The  total  H  field  in  the  plus  z  direction  is: 


Since : 


we  see 


H  =  (1  -  0)  k  -  0  M 


M 


-  H 


H 


0  _B_ 

1  -  0 


k 
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This  leads  to  the  load  line  construction  shown  in  Figure  A -8 


3-63-3864 


Figure  A -8  Load  Line  Construction 


4.  Relation  to  Cylindric  Rod 

Figure  A -9  shows  a  cylindric  rod  and  an  approximating  ellipsoid. 
If  we  select  the  dimensions  of  the  ellipsoid  so  that: 


b 

r 


Tfl 


then  the  two  volumes  will  be  equal. 


1.114 


3-63-3865 


Figure  A-9:  Fitting  an  Ellipsoid  to  a  Cylinder 
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The  fit  is  pretty  good  *  the  ellipsoid  protrudes  only  about  11  per  cent  out¬ 
side  the  cylinder. 

To  approximate  to  the  factor  f  employed  in  the  body  of  the  paper  we 
equate  the  slopes  of  the  two  load  lines  in  Figures  l  and  A -8 


f 


f 


2 

7T  r 


e  #  4 

1-8  n 


For  large  Awe  may  approximate: 


Log  Z  A  -  1 

A* 


so: 


4 

7T 


Log  Z  A  -  1 
,  Log  ZA 

A2 


A  typical  rod  with  a  length  of  50  cm  and  diameter  of  .  5  cm  has  a  A  of  100 
and  an  f  of  approximately  5.  5. 
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APPENDIX  B 


Air  Core  Loop  Relations 


The  following  formulae  are  useful  for  comparison  purposes.  In  all 
cases  we  have  been  concerned  with  the  loop  has  consumed  enormously  more 
power  than  the  rod  with  a  ferromagnetic  core.  Nonetheless,  the  formulae 
below  are  provided  to  facilitate  comparisons  in  cases  with  very  different 
limitations . 


NiA  =  7T  Ni  r" 


where  r  is  the  radius  of  the  loop.  Since  the  mean  turn  length  is  2  tt  r,  the 
wire  area  is  WV/ 2.71  Nr.  Hence  the  resistance  is: 


R 


coil 


WV 


Nr)2 


From  these  one  derives  the  hold  power: 


P 


hold 


P 

WV 


> 


Clearly  one  should  use  as  large  an  r  as  possible.  Moreover,  if  weight  is 
fixed,  one  wants  a  material  with  a  low  pd  product  -  i.e.,  aluminum. 
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